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1  Introduction 

J  - 

The  MODES  model,  discussed  in  PDRC  Report  84-06, Consists  of  the  LIFT- 
CAP  and  MRMATE  models.  LIFTCAP  is  a  network  flow  with  side  con¬ 
straints  model  while  MRMATE  is  a  network  flow  model.  MODES  uses 
Benders’  decomposition  as  a  solution  technique  that  separates  MRMATE 
and  LIFTCAP  and  adds  constraints  at  each  iteration  to  the  LIFTCAP 
model.  Benders’  thus  represents  a  finitely  convergent  algorithm. 

This  report  presents  an  alternative  procedure  that  maintains  problem 
structure  of  both  LIFTCAP  and  MRMATE.  Problem  sizes  are  maintained 
at  each  iteration  of  the  solution  procedure.  It  ensures  asymptotic  con¬ 
vergence  to  a  global  optimal  solution.  Due  to  its  nature  it  is  termed  the 
^averaging  method”?  The  steps  involved  in  implementing  this  procedure 
are  discussed  and  it  is  compared  to  Benders’  decomposition  technique. 

Computational  results  are  also  presented.  These  reflect  the  performance 
of  averaging  versus  Benders’  decomposition  technique.  - 

2  An  averaging  algorithm  for  MODES 

The  averaging  algorithm  for  MODES  works  as  follows 

1.  Solve  LIFTCAP  with  zero  costs  on  arcs  (or  some  reasonable  esti¬ 
mates)  and  generate  a  LIFTCAP  set  of  channels  co- 

2.  Solve  MRMATE  with  capabilities  c0  to  generate  a  solution  xj  and  a 
set  of  duals  u^v1  on  the  channels  and  MRs  respectively.  Set  t  =  1. 

3.  Set  the  costs  =  u*  on  the  channels  in  LIFTCAP  (there  is  one  dual 
value  associated  with  each  channel).  Solve  LIFTCAP  to  generate  so¬ 
lution  C{.  A  lower  bound  to  the  problem  is  given  by  LBt  =  LIFTCAP 
objective  value  -I-  £  MR  *  v*. 


='  =  £7  =  0 


>"'  +  T 


Solve  MRMATE  with  these  capabilities  c*.  The  objective  function 
value  provides  an  upper  UBt  bound  to  the  problem.  The  solution  to 
MRMATE  yields  duals  ut,vt. 


ilXU*  If  UM  Wl  WUWJV  mwwwtjwjw  w  u* wvtwuw!  w  yn  'jy  :jv.  'j*v*]F^s*xurr^runjr*jr*wwnijrw 


5.  Set 


<‘‘  =  ±ui  =  (i-]w-'  +  -t 


ut 


and 


‘  v<  .  1, 


vt 


6.  If  UBt  —  LBt  <  epsilon,  then  stop  with  channel  capabilities  cl  and 
allocations  xt.  Else,  go  to  step  3. 


While  implementing  the  algorithm,  store  the  maximum  among  lower 
bounds  generated  and  the  minimum  among  upper  bounds  generated  across 
iterations. 


2.1  Justification  for  the  averaging  procedure 


The  averaging  procedure  presented  in  Section  2.1  is  essentially  that  of  fic¬ 
titious  play  of  a  two  person  game.  Its  asymptotic  convergence  relies  on 
the  proof  of  [Robinson]  for  the  two  person  game.  This  section  discusses 
how  the  LIFTCAP  and  MRMATE  interaction  can  be  considered  as  a  two 
person  game.  LIFTCAP  attempts  to  generate  a  set  of  channel  capabil¬ 
ities  that,  when  passed  to  MRMATE,  yields  the  minimum  penalty  MR 
allocation  to  channels.  However,  since  the  LIFTCAP  feasible  region  is  a 
polytope,  the  particular  channel  capabilities  required  can  be  expressed  as 
a  convex  combination  of  basic  solutions  to  LIFTCAP.  Now  consider  the 
dual  feasible  region  of  MRMATE  for  a  given  LIFTCAP  solution.  Since 
it  is  unaffected  by  channel  capabilities  generated  by  LIFTCAP,  any  dual 
solution  can  be  expressed  as  a  convex  combination  of  dual  extreme  points 
(assuming  a  bounded  dual  region  achieved  by  the  use  of  dummy  nodes  in 
MRMATE). 

Let 


U  =  dual  feasible  region  of  MRMATE, 

C  =  primal  feasible  region  of  LIFTCAP, 
u*  =  dual  extreme  points  of  MRMATE, 
cl  —  primal  extreme  points  of  LIFTCAP, 


fk  and  gi  =  the  multipliers  associated  with  extreme  points  uk 
and  cl. 

The  objective  function  value  of  MRMATE  is  max^y  uc 

— ~  nisLXy^^  minc(c  tic 
=  max/  miny  E  Eg  fkUk  *  gic‘ 

=  max  min  E  E  fkgiukcf 

The  decomposition  of  LIFTCAP  and  MRMATE  can  be  considered  as 
a  two  person  game.  Given  all  its  extreme  points  u*,  LIFTCAP  generates 
a  mixed  strategy  g  that  minimizes  its  objective  function  value.  Similarly, 
given  extreme  points  c*,  MRMATE  generates  a  mixed  strategy  f  that  maxi¬ 
mizes  the  dual  objective  function  value.  This  produces  a  two  person  mixed 
strategy  game  with  extreme  points  of  the  dual  of  MRMATE  being  the  op¬ 
tions  of  player  1  and  extreme  points  of  the  primal  of  LIFTCAP  reflecting 
the  options  of  player  2. 

[Robinson]  has  shown  that  fictitious  play  of  a  two  person  game  con¬ 
verges  asymptotically  to  the  overall  optimal  solution  to  the  problem.  The 
algorithm  in  Section  2  is  precisely  fictitious  play  of  a  game  applied  to  the 
LIFTCAP/MRMATE  game.  Therefore,  the  averaging  algorithm  converges 
to  the  overall  optimal  solution  of  MODES. 

3  Comparison  with  Bender’s  decomposition 

3.1  Problem  Size  and  Structure 

The  averaging  procedure  maintains  problem  sizes  for  both  MRMATE  and 
LIFTCAP.  Thus,  MODES  could  be  run  for  as  many  iterations  as  required 
without  problems  of  memory  availability.  This  is  in  contrast  with  Ben¬ 
ders’  decomposition  which  adds  constraints  to  LIFTCAP  thereby  increasing 
problem  size  at  each  iteration.  Benders’,  in  current  implementations,  can 
be  run  for  only  about  10  to  20  iterations  before  numerical  stability  prob¬ 
lems  prevent  furthur  iterations  with  the  network  flow  with  side  constraints 
LIFTCAP  solution  code. 

Furthermore,  the  LIFTCAP  model  structure  is  maintained  across  iter¬ 
ations.  Channel  capabilities  generated  at  each  iteration  have  an  intuitive 


s  s. * * 
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interpretation  as  the  best  possible  asset  assignment  for  the  given  set  of  costs 
(MRMATE  dual)  allocated. 


Flexibility 


The  averaging  algorithm  works  for  any  set  of  weights  (0  <  /  <  1)  used  to 
create  a  new  set  of  channel  capabilities  or  MRMATE  duals  at  each  iteration. 
Thus,  the  sequence  of  channels  selected  can  be  controlled  by  adjusting  the 
weighting  parameter. 

If  MODES  is  initialed  with  a  good  set  of  asset  allocations,  then  the 
weights  on  the  new  sets  of  channels  could  be  set  very  small.  As  the  iterations 
progress,  however,  the  weighting  parameter  could  be  increased  to  introduce 
better  channel  assignments. 


3.3  Computational  Results 


An  example  problem  was  used  to  produce  test  results  from  the  use  of  Ben¬ 
ders’  and  averaging.  The  problem  consists  of  100  nodes  in  LIFTCAP  with 
1600  variables  and  an  MRMATE  problem  that  has  150  MRs  and  300  chan¬ 
nels.  Figure  1  shows  upper  bounds  generated  at  each  iteration  by  the  two 
procedures. 

Figure  2  shows  lower  bounds  generated  by  averaging  and  Figure  3  shows 
upper  bounds,  while  Figure  4  shows  upper  and  lower  bounds  for  averaging 
across  a  hundred  iterations. 

Figure  5  gives  the  fraction  of  MRs  allocated  to  channels  for  the  two 
procedures.  The  averaging  procedure  appears  to  perform  well,  probably 
because  it  results  in  more  non  zero  channels  very  quickly,  thereby  improving 
the  MR  allocations  as  compared  to  Benders’.  This  is  another  desirable 
feature  of  this  algorithm.  Continued  testing  of  this  procedure  is  necessary 
before  strong  conclusions  can  be  made  regarding  averaging  versus  Benders’. 


4  Conclusions  and  extensions 

/ ,  ..  7 

This  report  outlines  a  procedure  to  solve  MODES  without  adding  Benders’ 
constraints  to  the  LIFTCAP  model.  -The  procedure  is  asymptotically  op¬ 
timal  and  maintains  problem  size  and  structure  at  each  iteration.  Some 


!,%»*  ,»**.»**  VI,  Vf. 


U  »*»»**.  >*!.«*>  gMgl  **» ,tT.|M 


h'M 


Brag 


m. 

"V 


m 

Bjjft 

iKt't 


Vm 


i 


300 


!*■*»•*•»*  »»b  |.|  »,%  i.l 


I  4  |  .u  ViyiVtWJjVt^a  *  fa*  L*  ’  fi*  left  J»*  JtMjV|V|Vl  Vl 


prelixnnary  results  were  presented  which  demonstrate  the  effect  of  this  pro¬ 
cedure  over  Benders’  decomposition  technique.  '  -  -  „  - 

Further  research  is  in  progress  to  determine  suitable  choices  of  weighting 
parameters  at  each  iteration.  Furthermore,  the  most  suitable  procedure  for 
MODES  might  be  a  hybrid  approach  that  combines  a  few  iterations  of 
Benders’  with  the  averaging  procedure. 
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